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Abstract 

The S-meson decay constant can be measured on the lattice using a 1 / m& ex- 
pansion. To relate the physical quantity to Monte Carlo data one has to know the 
renormalization coefficient, Z, between the lattice operators and their continuum 
counterparts. We come back to this computation to resolve discrepancies found in 
previous calculations. We define and discuss in detail the renormalization procedure 
that allows the (perturbative) computation of Z. Comparing the one-loop calcu- 
lations in the effective Lagrangian approach with the direct two-loop calculation 
of the two-point S-meson correlator in the limit of large 6-quark mass, we prove 
that the two schemes give consistent results to order a s . We show that there is, 
however, a renormalization prescription ambiguity that can have sizeable numer- 
ical consequences. This ambiguity can be resolved in the framework of an 0(a) 
improved calculation, and we describe the correct prescription in that case. Finally 
we give the numerical values of Z that correspond to the different types of lattice 
approximations discussed in the paper. 
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1. Introduction 



The study of heavy flavours has attracted a lot of interest in the last few 
years, both from the experimental and theoretical point of view. Many interesting 
properties of matrix elements can be derived on rather general assumptions, using 
the spin symmetry that emerges in the limit in which the mass of the heavy quark 
goes to infinity [1], [2]. 

Beyond symmetry arguments, the only method based on first principles that 
can be used to predict physical quantities, such as the S-meson decay constant, 
or the universal form factor for current matrix elements between heavy mesons, is 
lattice QCD, complemented with a l/m^ expansion of the heavy-quark propagator, 
as proposed in ref. [1] . 

In this scheme, the heavy-quark field is no longer explicitly present on the 
lattice, but it is replaced by a static colour source, and one has only to consider 
correlation functions with gauge and light-quark fields. The latter ones can be 
evaluated in lattice QCD by Monte Carlo simulations in the standard way. This 
method is the only one applicable to the B meson since the b quark is too heavy 
to live even on lattices with the smallest lattice spacing presently attainable in 
numerical simulations (a -1 of the order of 3 to 4 GeV). In the phenomenological 
analysis, one then has to interpolate between the infinite mass limit computed in 
this way and the results obtained with the standard approach for quarks lighter 
than a" 1 [3]. 

From now on we will concentrate on the method based on the infinite mass 
limit. We will call it for short in the following the static quark approach. 

As usual, we have to deal with the problem of lattice renormalization to relate 
numbers effectively measured on the lattice to the corresponding physical quantities. 
Two methods have been used for the calculation of the necessary renormalization 
constants [4], [5]. Apparently the two approaches have led to different numeri- 
cal results. One of the purposes of this paper is to understand and to solve this 
discrepancy. 

Let us focus on the computation of f B , the decay constant of a S-meson. From 
the current-current correlation function: 
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we can extract the value of f B by looking at its large t behaviour: 
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where Mb is the B- meson mass. 

The presence of the heavy b field prevents a direct lattice evaluation of the 
expectation value in (1.1). The nice idea of ref. [1] is to replace, in the limit — > oo, 
the b and b fields in eq. (1.1) by the static b propagator: 



S 6 (x', 0; x, t) = V x l^j 8(x - x') [#(£) e 
where m& is the 6-quark mass 1 and 



, (1-4) 



= Pexp 
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is a path ordered product along a temporal line with T a the usual Hermitian colour 
matrices. One formally gets in this way (t > 0) : 



P(f) 



P stat {t) = e~ m ^ (Tr 



^~Y^ 7075 9(0,t) 9(0,0) 7075 ^oQ}> • (1-6) 



Except for the c- number factor e~ mbt , the operator whose vacuum expectation 
value is taken in (1.6) depends only on light quarks and gauge fields. It is a gauge- 
invariant quantity, as expected. On the lattice the Green function corresponding to 
the bracket in eq. (1.6) is: 

P L (r) = ( Tr(^^Sl(0, r; 0, 0)U (l)U (2)...U (r - 1))) (1.7) 

where r = t/a, U^{x) is the link matrix pointing from the site x in the direction fi, 
and Si is the light-quark lattice propagator. The naive a — > limit transforms (1.7) 
into (1.6) and thus a Monte Carlo estimate of (1.7) would give an estimate of (1.1), 
i.e. of f B . However, the breakdown of the static approximation for the propagator 
(1.4) at high frequencies (larger than m&) and, as usual, the lattice granularity 
induce renormalization effects, so that we expect a relation of the type: 



P{t) 



a.T = t — * oo 



ZPmc(t) e~ EBa \ 



(1.8) 



1 As we shall see, we do not need a precise definition of 
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where Pmc is the numerical value of Pl [eq. (1.7)] as measured in Monte Carlo 
simulations; Eb is introduced to match the time dependence of the two sides of (1.8); 
its role and significance will be discussed in more detail later. The renormalization 
constant Z is computed by comparing in perturbation theory the corresponding 
continuum and lattice quantities 2 . 

To this end two methods have been proposed. The first one (BLP method) [4] 
consists in a direct perturbative comparison of the correlation functions (1.1) and 
(1.7). To order g 2 this amounts to evaluating the two-loop diagrams of Fig. I 3 . It 
will also be convenient to compute (1.6) as an intermediate step, which involves 
considering the graphs of Fig. 1 in the case of a static continuum quark source. The 
contributions to Z coming from the diagrams of Fig. 1 labelled by 6, c and d will 
be called Zu g ht, Z vertex and Z heavy respectively. To the order we work we have 

Z = ZlightZyertexZheavy ■ (1-9) 

The EH method of ref. [5] relies on the remark that the dynamics of a 
static (colour triplet) quark source can be represented by the effective Euclidian 
Lagrangian: 

C E h = ~& (do + igA )B , (1.10) 

where the B (B^) field annihilates (creates) a static quark. It should be noted 
that a fixed four- momentum (m&, 0) has been removed from the momentum of the 
heavy quark. A completely independent static field should be introduced to describe 
processes involving heavy antiquarks. The B field is a two-component spinor that 
describes spin-up and spin-down static b quarks. 

The free propagator in position space is given by: 

S EH (x,y) = (B(x)B\y)) = e(x -y ) <5(x - y) , (1.11) 

2 For simplicity, we will use a massless light quark, since none of our conclusions 
depends on this choice [4], [5]. The Feynman gauge is used everywhere. Con- 
tinuum renormalizations have been performed, as in ref. [4], in the MS scheme 
supplemented by an on-shell mass renormalization. 

3 In order to avoid redundancy, we use the same graphic representation for the 
three cases we shall encounter. According to the situation, the heavy line must be 
interpreted as a continuum quark propagator (1.1), a lattice static quark source 
(1.7) or a continuum static quark source (1.6). The context should make it clear 
which interpretation is to be used in each case. 
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and in momentum space by: 



Seh( P ) = — ! — ■ (1.12) 
Po + it 

Equation (1.10) can be formally obtained starting from the Dirac action of a particle 
of mass nib and integrating out the "small components" of the Dirac field in the 
limit nib —> oo. 

The lattice-discretized version of the Lagrangian (1.10) can be chosen in Eu- 
clidean space to be [5]: 

C L EH = B\n) (B(n)- [t/ (n-6)] t B(n-0)) , (1.13) 

with n = (n, r), leading to the free propagator for the lattice static source: 

s L EH {p)= i( — - 7T-. ■ (1-14) 

In an external gauge field the exact static-quark propagator coming from the action 
(1.10) [resp. (1.13)] leads indeed to eq. (1.6) [resp. (1.7)]. 

In the presence of a "local mass" term MaB^(n)B(n) the free propagator com- 
ing from (1.13) becomes 

<B(n)Bt(0) >=e- Ma ( T+1 ^(r)5 n , . (1.15) 

The main reason to introduce an effective action for the heavy quark is to reduce 
the determination of the renormalization coefficient Z to the usual calculation of 
vertex and self-energy corrections, i.e. of the diagrams of Fig. 2. In this way, to 
order g 2 , it will be sufficient to compute one-loop graphs, whereas, in the BLP 
method, one has to deal with two-loop calculations. 

It should be noted however that one has to face, in the effective-action method, 
the problem of principles that the correlation function we are interested in appears 
as a product of two operators at different times, but at the same space position. 
Since the effective theory is not Lorentz-covariant, the renormalization procedure 
is not straightforward and it is not obvious that complications, coming from an 
interplay between infrared singular behaviour in time and ultraviolet divergences 
in space, will not be present. In this work we shall prove the important fact that 
indeed such complications do not arise in the one-loop calculation of Z. In spite 
of that, the two papers [5] and [4] give different results. The second purpose of 
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this paper is then to study in detail the nature of the normalization condition that 
determines Z. 

We will see that, as a matter of fact, the difference between the values of Z 
quoted by different authors is due to terms that originate, in the language of ref. [5] , 
from the freedom in the choice of the expression for the effective heavy-quark action 
on the lattice, while in the language of ref. [4] they can be seen to come out from an 
intrinsic (multiplicative) ambiguity in the form of the normalization condition that 
determines Z (see also [6]). Although these terms vanish as a — > 0, in the present- 
day simulations they are numerically not negligible. We shall prove, however, that 
problems related to this ambiguity can be eliminated, to leading order in a, by a 
consistent use of an 0(a)-improved action for the light quarks [7], together with the 
appropriate form of the normalization condition. 

The paper is organized as follows. In section 2 we discuss in detail the nature 
of the normalization condition that allows us to relate, for large t, the physical 
Green function P(t) to the Monte Carlo measured quantity Pmc(t)- The detailed 
comparison of the approaches of refs. [4] and [5] is presented in section 3, where the 
controversial calculation of the self-energy of the heavy-quark line is discussed. We 
show that numerically the difference between the two approaches lies in the way the 
"wave-function" renormalization constant of the heavy line is extracted from the 
relevant Feynman diagrams. In particular we notice that in the scheme of ref. [5] the 
result depends upon the way the heavy-quark action is discretized on the lattice, 
section 4 addresses the problem of the potentially dangerous interplay between 
infrared (in time) and ultraviolet (in space) singularities in the vertex diagram. We 
rigorously prove that, to the order at which we work, this phenomenon does not 
occur, and the two methods lead to the same result for Z vertex . Along the way we 
correct an error present in ref. [4]. In section 5 we give the final numerical values 
of Z. 



2. The normalization condition 

The value of Z in (1.8) can be obtained from the comparison of the large t 
behaviour of the non-perturbative ratio 

<J>,0) J qb (x,t)>\ phya Z&Le-Mst 

K NP — „ ~~Ft -m 

Pmc{t) » t = 1 ^ 00 OMce 



•5 



to the corresponding perturbative one 



< 4,(0,0) J qb (^t)>\ PT 



p L (t)\ 



PT 



In eq. (2.1) B is a mass parameter fitted from the (expected) time exponential 
behaviour of Pmc( t ) : 

Pmc(t) C MC e~ m ■ (2.3) 



To any order in perturbation theory the behaviour of the r.h.s. of eq. (2.2) is, 
actually, a polynomial in t. However in the large t limit, this polynomial is the g 2 
expansion of a single exponential. In the following, for brevity, we will simply refer 
to it as an exponential. 

With this in mind, one can evaluate Z through the following two-step procedure: 

1) given the quantity, indicated as Rnp in (2.1), one computes in perturbation 
theory exactly the corresponding ratio of Green functions, Rpr] 

2) once the large t behaviours of Rnp and Rpt have been removed, the two 
remaining constant quantities are declared to be equal in the limit of vanishing 
lattice spacing. In formula (2.2), having defined Mpt so as to perturbatively 
remove the exponential t behaviour of Rpt, i.e. 

lim e MpTt R PT (t) = const , (2.4) 

the normalization condition can be expressed as: 



— ^— = lim lim e Rpx{t) = Z . (2.5) 
For future reference we write: 

M PT = dpT-a , (2.6) 

where, by definition, a is the perturbative part of B: 

B = a + S NP . (2.7) 

We note that, although a is, as we shall see, linearly divergent as a — > 0, the 
non-perturbative quantity Snp is regular, i.e. 

lim B — a= lim 5mp = finite (2.8) 

V ' 



Two observations are in order here. The first one is that the normalization 
condition (2.5) has to be imposed in the continuum limit, i.e. among finite quantities 
in the limit a — > and ar = t fixed. The second, and more important one, is that 
eqs. (2.4) and (2.5) can fix Z only up to a multiplicative factor (going to 1 as a — > 0). 
In fact the condition (2.5) does not forbid multiplying Rnp by a time-independent 
constant, provided we similarly multiply Rpr by the corresponding perturbative 
factor, as prescribed by the normalization procedure described above. For instance 
we can equally well replace eq. (2.5) by: 

f l MB p M PT t p l,\ 

lim 1^ = li m lim £ tlElStl (2.9) 

a^O e~ xaB CMC ar=t^oo g-xacr ^ > 

where x is any (small) real number. Equation (2.9), for x = 1, is essentially the 
normalization condition imposed in ref. [5] (see section 3), while for x = we get 
the prescription of ref. [4]. 

The exponential exp(— xaa) in the r.h.s. of eq. (2.9) is introduced to com- 
pensate for the corresponding factor, exp(— xaB), present in the l.h.s. This is in 
agreement with our prescription and it is necessary because, as we have said, B is 
(perturbatively) linearly divergent, so aB does not vanish as a — > 0, while a(B — a) 
does. 

Using eq. (2.5) or (2.9), one can extract f B from Monte Carlo data. The reason 
why in practice the two choices, x = and x = 1, lead to different values for f B 
is then obvious. Since Pmc and, consequently, B is only known at finite values of 
a, the two exponentials in eq. (2.9) do not compensate exactly. If a is taken to 
one-loop in perturbation theory, one has: 

exp[xa(B-a)] = l + x[0(a)+0(g 4 )] ^ 1 . (2.10) 

As a consequence, the actual value of f B , extracted from Monte Carlo data, depends 
on the detailed form of the normalization condition. From eq. (2.9) one gets: 

(f B )l = (f B )l =Q eMMB-cr)] (2.11) 

where, using the definition of Z given by eq. (2.5), one has 

^ (f B )Lo = ZC MC • (2.12) 
For x = 1, eq. (2.11) becomes 

(/b)e« = B ^ BLP eX P ' 



(2.13) 
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which gives the relation between the value of f B extracted from Monte Carlo data 
following the approach of ref. [5], called here (f B ) EH , and (/s) BIiP , the same quan- 
tity as evaluated within the approach of ref. [4]. 

Strictly speaking, there is no compelling reason to prefer x = to x = 1, or, for 
that matter, to any other value of x. We shall see, in fact, in subsection 3.4 that 
exactly the same line of arguments, which "naturally" led to the choice x = 1 in 
ref. [5], can also be used to arrive "naturally" at the value x = — 1. 

The question is then whether there is any criterion to fix this problem univo- 
cally. It is not difficult to see that this can be done if we decide to improve the 
light-quark Wilson action [7] employed in the Monte Carlo simulation of Pmc an d 
in the corresponding perturbative calculations of Z [8], [9]. 

In this case, in fact, since eq. (2.8) maintains its validity, the only consistent 
way of eliminating all 0(a) terms from the lattice evaluation of f B is to take x = 0, 
that is to use the normalization condition (2.5). If we do so, we remain with only 
0(ag 2 ) + 0(g 4 ) corrections. 

3. Heavy-quark renormalization 

In this section we want to compare in detail the approaches of refs. [4] and [5] 
and to clarify the origin of the discrepancy between the two results by discussing the 
calculation of the graphs of Figs. 1 (d) and 2 (c) which give rise to the controversial 
contribution, ^ eauy , to Z. 



3.1. BLP evaluation of Zh 



eavy 



With the BLP method the graphs of Fig. 1 (d) have been computed in ref. [4]. 
In the Feynman gauge, their contributions to (1.7), (1.6) and (1.1) are respectively: 



Ptd) - / - 



d 4 k l — cos(kot/a) 



ytree 
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(3.1) 



- aar \ P tree (t) 



P&atV) = ± 



3tt 2 
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d( + \ _ „-m b t _9_ 



P d (t) = e 
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ytree 



(t) 



(3.2) 
(3.3) 



where, on the lattice: 



Tjtree 



V d 4 k 1 
{t) = J (2^ 2 



x Tr 

ptree 



(l + 7o) 
(f) . 



E M 7^ sin p M + ma + r - cos p u ) 
E M sin p M + (ma + rJ2u i 1 ~ cos p v )) 



(3.4) 



and in the continuum 



jtree 



In (3.1), -Foooo and -Foooi are the numerical constants defined in ref. [10]: 
Fqooo — 4.369 and Fqooi — 1.311; = 0.5772 is the Euler constant, and 



J2 - +7r 



a J 



24tt 3 a J E(l - cos fc,-) , 

t ' (3.5) 

• 6.65 

47T 2 a 

is the linearly-divergent one-loop radiative correction to the energy of the static 
source. In the following we will often write 1 — aa = e~ aa + 0(g 4 ). 

From these results we get for the Zheavy contribution to the total normalization 
constant Z, defined in (1.9) and (2.5): 

Z heavy - 1 = lim e M ^ PtrCe{t) + Pd(t \ - - 1 



Pt^(t) + Pi{t) 



9 L ( u l, (m b \ is F 0000 + ^oooi J (3.6) 
g2 ^ 4 ln(m 6 a) + -lnf^^-2.84) , 



4tt 2 [3 v ' 3 \ /I 

where e Mdt is the appropriate factor necessary to compensate for the overall time 
exponential behaviour in the r.h.s. of the first line of eq. (3.6). In perturbation 
theory the multiplication by this factor simply amounts to subtracting out from the 
computed amplitudes all terms proportional to t. 
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3.2. EH evaluation of Zheavy 

With the second method [5], the self-energy of the heavy quark, T,(po) EH , is 
computed directly with the Feynman rules deduced from the effective action (1.13), 
obtaining: 

4 „i + r n d 4 k 



2 ± / A(k +2ap ), 



3 a_7 (2tt) 4 4£sin 2 (%) e J ( fc °+ a Po) - 1 - e 

(3.7) 



2 1 V r/ 4 k i 



3 y a _7 (27r)4 4Esin 2 ( ^ ) 



2 



It is crucial to remark that if we insert T,(p ) EH into the external quark loop, 
we recover exactly the integral giving in eq. (3.1). So the two methods should 
agree. In ref. [5] , the wave function renormalization was computed with the standard 
formula: 

It is immediate to check that the radiative correction to the mass is the same as in 
(3.5), i.e. 

a = £(p = 0) . (3.9) 

Using eq. (3.8), the result for the ratio of the wave function renormalization 
constant of the continuum dynamical 6-quark [eq. (1.1)] to that of the lattice static 
B field in (1.13) is: 



dpo 



dp 

with Z heavy given by eq. (3.6) 



= Z heavy (1 - aa) ~ Zheavy e aa , (3.10) 

Po=0 



3.3. Comparison of the two approaches 

To compare the results for Z obtained in refs. [4] and [5], one has to include 
in the calculation the values of Zu g ht and Z vertex . As for Zu g ht, it is easy to see 
that the analog of the ratios in eqs. (3.6) and (3.10) for the light-quark self-energy 
diagrams, explicitly given in subsection 5.1 [Figs. 1 (b) and 2 (a)], yield the same 
results, both with the method of ref. [5] and with that of ref. [4] . The reason is that 
there is no ambiguity of the type we have discussed above for the heavy b quark. In 
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fact, on the one hand the way the mass-counterterm for the light quark is introduced 
is univocally fixed by the form of the Wilson action that is being used, and, on the 
other hand, both in the Monte Carlo simulations and in the perturbative calculations 
the same expression for the light-quark action with the same mass renormalization 
condition (vanishing pion and quark masses respectively) is employed. Concerning 
the contribution of the vertex diagrams [Figs. 1 (c) and 2 (b)], we will show in 
section 4 that the results of refs. [4] and [5] for Z vertex also coincide. 

It then follows that the relation between Z in (1.8), as computed in ref. [4] and 
the current renormalization constant Zj introduced in [5] is 

Zj 2 = Z e~ ao ". (3.11) 

Now we turn to the central issue: given a Monte Carlo data set, what is the pre- 
diction for the physical f B , and do the two methods agree? 

Fitting the Monte Carlo data for large t through formula (2.3), in ref. [4] one 
arrives at the result 

B ' = ZCmc ■ (3.12) 



\BLP 



In ref. [5], the Monte Carlo data are fitted instead, in view of eq. (1.15), ac- 
cording to 

P MC ~ A e ~ B ^ . (3.13) 

t — >co 

Comparing (2.3) and (3.13) one has 

C MC = Ae~ Ba . (3.14) 

For the physical correlation function (eq. (26) in [5]) one obtains: 



f -^ ± \ EH e~ MsaT = Zj 2 e- mbaT A e -Ba(r+l) e -Aa(r+l) _ (3^) 

2 



For consistency of our notation we have rewritten the factor 1/(1 + a5m/Z) of ref. [5] 
as e~ aA . By matching the large t behaviour of (3.15), one gets: 

A = M B -m b -B (3.16) 

and 2 

f -^-\ EH =ZJ 2 Ae- Ba e- Aa = Zj 2 Ae- {MB - mb)a . (3.17) 

At this point the authors of ref. [5] identify Mb — mi, as the S-meson binding 
energy. They estimate A in the limit a — > 0, by arguing that Mb — must 
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converge to a constant in the continuum limit, i.e. lima(Ms — m&) = 0, so they end 
up with: 

£^l.„ = A . (3.18) 

Remembering the relations (3.11) and (3.14), the comparison of (3.12) and 
(3.18) gives: 

a(B-q-) 

JbIeh = JbIblp xe 2 • (3.19) 

This answers our central question: given a set of Monte Carlo data the number 
extracted for the physical f B according to [5] [eq. (3.18)] will be that of [4] [eq. (3.12)] 
multiplied by exp[ a ^ jB 2 ~ (J - > ]. 

Since a(B — a) = a(— A — a + Mb — rrib) is 0(g 4 ) + O(ag ), it follows that both 
results agree in principle in the approximation in which we are working. However, 
taking a typical Monte Carlo slope aB = 0.7 (aB = 0.5) for j3 = 6 (/? = 6.4), eq. (3.18) 
gives a value for f B \ EH that is larger by a factor ~ 1.3 (~ 1.2) than f B \ BLP for the 
same Monte Carlo data set: the discrepancy has an obvious practical relevance and 
its origin should be fully understood. 



3.4. Further remarks 

Let us insist once more on the content and the origin of the factor in (3.19). 
First, it has nothing to do with the computational method, namely the use of either 
direct computation of the correlation function or effective Lagrangians. Second, 
there is no discussion on the validity of the starting equation (1.8). Thirdly there 
is no discrepancy in the mathematics. 

The point is that, in ref. [4], Z is computed (to second order of perturba- 
tion theory) by using the normalization condition (2.5), thus directly obtaining 
eq. (3.12). In this way, one must simply assume that an Eb exists, which matches 
the time dependence of both sides of eq. (1.8), without the drawback of having to 
give a problematic estimate of a linearly divergent quantity. As we have discussed 
in section 2, however, an apparently irrelevant modification of the normalization 
condition, such as the one made in eq. (2.9), leads to a non-negligible numerical 
ambiguity in the estimated value of f B given by the factor (2.10). 

The approach of ref. [5], instead, relies on the use of an effective Lagrangian for 
the heavy quark, both in the way the large t dependence of the perturbative and the 
Monte Carlo expressions of Pl is parametrized and in the way the wave-function 
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renormalization of the heavy-quark line is computed. The authors of ref. [5] were, 
in fact, led by the form (1.15) of their free propagator to the parametrization (3.13); 
they were also driven by formula (3.8), which looks so familiar, to choose Z 2 as the 
quantity to be computed perturbatively. 

The key observation at this point is that the effective action is merely an 
intermediate step used to compute the relative normalization between the physical 
observable and the quantity measured on the lattice. Neither quantity has anything 
to do with any effective action and the final result should not depend on it. However, 
it is easy to produce a variety of effective lattice actions whose continuum limit is 
(1.10) and which, by the use of the corresponding matching condition and of formula 
(3.8), lead to different values of f B for a given set of Monte Carlo data. 

For instance instead of (1.13), one could as well choose for the discretized 
version of (1.10) the Lagrangian: 

C' EH = B f (n) (t/ (n + 6)B(n + 6) - B(n)) . (3.20) 

If a local mass term MaB> (n)B(n) is added, one gets: 

<B(n)B t (0)>=e- Ma(r - 1 ^(r)5 n , o . (3.21) 

Notice the change (r + 1) — > (r — 1) in eq. (3.21) with respect to eq. (1.15). 
This change propagates step by step from eq. (3.13), leading now to: 

f' B \ EH = / B LpXe"^ (3.22) 

instead of (3.19). 

Another interesting choice is to use an (9(a)-improved effective Lagrangian for 
the heavy quark, i.e. to take: 

C imp = Bt(n) ( ^B(n) - 2£/ (n-6) t B(n-0) 

V , (3.23) 

+ -C/b(n-6) t C/o(n-26) t B(n-26)j . 

One checks that, if the static field is integrated out, the Lagrangian (3.23) leads 
exactly to the same correlator (1.7) as the Lagrangians (1.13) or (3.20), up to 
exponentially small terms, vanishing as t — > oo. 

In view of this argument, the authors of [8] and [9] claim correctly that there 
is no need to use an improved effective action for the heavy quark, if one uses 
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the normalization condition (2.5). However, if one insists on using eq. (3.8) to 
define Z heavy the precise form of the effective action does matter. In fact, the free 
propagator obtained from (3.23) is: 

S imp (p) = _i {2eipoa _ i ei2poa _ 3 } + - e • (3-24) 

It coincides with the continuum one [eq. (1.12)] up to 0(a 2 ), while (1.14) differs 
already at O(a). We have computed T, imp (p), the self-energy of the heavy line, with 
the "improved" effective Lagrangian (3.23) and derived the mass and wave function 
renormalization, by using formulae (3.8) and (3.9); £j mp (0) is identical to (3.5), but 
the noteworthy fact is that, instead of eq. (3.10), one obtains 



1 -i 



8T, C 



dpo 



.BY,. 



dpo 



'heavy 



(3.25) 



po=o 



Following the steps of ref. [5], we end up this time with a result identical to 
(3.12) (the same as in ref. [4]), that is: 

(/fllimp \eh (/b)x=0 /bIblp 

Thus a determination of Z that uses formula (3.8) can give the same result as (2.5), 
if one uses the Lagrangian (3.23). 

From all this analysis, we can conclude that the use of the effective action 
approach by itself is not sufficient to lead to an unambiguous (in the sense we are 
discussing here) determination of f B from Monte Carlo simulations. 

We have instead shown in section 2 that this can be achieved to 0(a) only if 
one consistently uses: 

i) an O(a) improved action for the light quarks. 

ii) the normalization condition (2.5), in which the matching does not involve un- 
known O(a) terms. 



(3.26) 



4. Vertex renormalization 

In this section we discuss the evaluation of the contribution to Z of the vertex 
diagrams [Fig. 1 (c)] and we rigorously prove that the effective Lagrangian method 
leads to the same result as the full-fledged two-loop calculation of ref. [4]. That is 
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to say, there is no dangerous interplay between IR and UV divergences. The direct 
evaluation [4] of the two-loop graph of Fig. 1 (c) gives: 

P[(t) = A -g 2 ^ (J[ att (r) + X l 2 att (r)) , (4.1) 

with 



+ 7T 



rlatt 



dpo 
2n 

-ipo 



+7TT 



d 3 ] 



— ITT 
+ 7T 



r cr\ 
1 72tt 



(2tt)* 



d 4 k 



1 ~e ipo + £(*) 
2A 2 (p , *) 



x 



Jk /2 



(2tt) 4 e lk ° -lAi(jfc) 



-ifco/2 



1 _ e i(po+fco) + £(E +k) 

A 2 ( Po + k , E + k) 



o/2 l-e l ^°- fc °) + S(E-k) 



A 2 (p -A;o,^-k) 



2*»« = 3£ T ^ 3 P e i P0 r -^in(^) x 



2tt 7 (2tt)3 2A 2 (;po,^) 

— 7TT 



+ 7T 



_, po / _^k_ e ifco/2 ! , „_, fcn/2 -isin^ + fcj) 



( 27r )4 e ifco - 1 Ai(fc) I A 2 (po + fc , * + k) 

ifco/2 



e 



/2 -isin^-fc,) 



A 2 (p -^o,^-k) 



(4.2) 



where the A's come from the gluon and fermion lattice propagators: 
Ai(fe) = ^(l-cos(fc A )) 

Z A 

A 2 (p ,p) = E sin2 (^) + ^ [1 - cos(pa)]^) • 



(4.3) 



As we said, we have chosen here to work with a massless light quark, since the 
renormalization constants are independent of the mass [4], [5]. No mass has been 
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given to the gluon since the quantities we are dealing with are infrared-finite after 
the singularities among individual terms have been cancelled. 

To determine Z vertex , we have to study the behaviour of T l £ u et J 2 a " for large t, 
at fixed a, i.e. when r goes to oo. The integrations over the temporal components, 
Po and ko, are performed exactly, by closing the contour in the complex plane. 
Oscillating factors are turned into exponentially damped ones, easier to handle. 
The large-r behaviour of the resulting spatial integrals is analysed by means of 
Lebesgue's lemma. We find (see appendix A for a sketch of the proof): 



T latt(-.\ _2_ f 

1 1 ' -co ^ J (27T 



1 " " ' ik ' 1 V, (1 - cos(jfej)) 1 - < 
(2tt) 4 2Ai(fc)A 2 (fc ,k) 1 - e ik ° 



+ d+o(ln(r)/r) , (4-4) 
l l 2 att (r) ~ C 2 + o(ln(r)/r) , 

T — »00 

where Ci and C 2 are some (r-independent) continuum-like integrals. A similar 
analysis can be performed in the continuum for -P s c tat (t) [eq. (1.6)], using dimensional 
regularization (d = 4 — e) to control its ultraviolet behaviour, and gives: 

P c statit) = \g 2 ^ + J5 to V)) , (4-5) 



with 



+oo 



d 4 ~ 6 k e ik(>t -l 

(2^PW ' W± ' " (4.6) 

J?°V) ~ C 2 +o(e) 



The constants Ci and C 2 are the same as in eq. (4.4). The decomposition into 
the terms appearing in eqs. (4.4) and (4.6) were performed in such a way that the 
dependence upon the regularization scheme was isolated in the first terms, leaving 
C\ and C 2 renormalization-scheme independent quantities. 

Using the MS renormalization scheme in the continuum, we find that the 
contribution to the ratio of (1.6) to (1.7) coming from this graph (see appendix B) 
is: 
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9 2 4 



- 1 



8tt 2 3 



+ 7T 

21„ W + 7E -F 0000 --/ rf »k s - I L^ 



+ 7T 



4tt 2 



/ £sin 4 % 



(4.7) 



d 4 k 



AfA 2 (£; ,k) 4A 2 (A; ,k) 



g 2 ( 4 ln(/m) 



4tt^ 



+ 7.79 



where the factor e Mct plays the same role as the similar exponential in eq. (3.6). 
The result (4.7) coincides with the value obtained from the evaluation of the one- 
loop vertex graph of ref. [5]. This indeed shows that there are at this order no 
complications due to the problem of the renormalizability of a non-covariant effective 
theory. In ref. [4] this computation was actually performed with a non-vanishing 
light-quark mass. The results agree, apart from the l/2ir J <i 3 kA 2 ~ 1 (0, k) term, 
which is present in eq. (4.7) but was erroneously forgotten in ref. [4]. 



5. Results 

In this section we complete our discussion on the evaluation of Z, by giving 
the expression of Zn g h t , and we compare the values of f B obtained with different 
normalization conditions, with the number one would obtain by using an O(a)- 
improved light-fermion action [in which case it is the form (2.5) of the normalization 
condition that, for consistency, has to be employed]. 



5.1. Contribution of Zi ight 

The last piece of the puzzle is the contribution to the renormalization constant 
due to the graph 1 (b) which is the usual wave-function renormalization for a Wilson 
fermion [11], [12], [13]. Numerically one obtains: 

where, once again, e Mbt is the appropriate factor introduced to make the ratio in 
the l.h.s. go to a constant as t — > oo. 
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To reach the final result, we need the contribution to the ratio P s tat/ Pl coming 
from the d-type graphs in Fig. 1 [in analogy with eqs. (4.7) and (5.1)]. One has: 



M P tree (t) + PLt(t) 

ptree {t) + pd {£) 



1 



g 2 f 4 Infra) , 2 
4^ + 3 7£ 



-^0000 + ^000 1 \ 

3 J 



9_ 

4tt 2 



Infra) - 1.51 J 



(5.2) 



The relation between the two continuum correlation functions (1.1) and (1.6) 
is already known [4], [5]; the result is the same in the two methods and it has also 
been checked in the temporal gauge ( A = 0) [14] . It can be written: 

=2/3" 



Mt 



P(t) 



stat 



(!) 



PT 



9 
2tt 2 



In 



lie 



(5.3) 



We insist once more on the fact that all the exponential factors in eqs. (3.6), 
(4.7), (5.1), (5.2) and (5.3) are defined so as to insure the constancy of the corre- 
sponding ratios of correlation functions as t — > oo. The values of the mass coefficients 
are totally irrelevant for the purpose of computing Z, because the only thing one 
actually has to do is to neglect in the calculations all the terms that are proportional 
to t. Combining eqs. (4.7), (5.1), (5.2) and (5.3), we obtain the result: 



M PT t 



P(t) 



P L (a,t) 



= Z 



1 = 



PT 



4tt 2 



21n( 



nibae 



-2/3 



13.59 



(5.4) 



5.2. Numerical evaluation of Z 

The relation of the physical value of f B to the measured lattice quantity 



rlatt _ 

J B 



[2Cmc \ 



1/2 



V M B ) 

extracted from the behaviour of Pmc(j) at large t [eq. (2.3)] is then, from eq. (5.4): 

,2 



fs=f 



latt I -i , 9_ 

B + 8tt 2 



2 In (. 



arribe 



-2/3 



13.59 



(5.5) 



Numerically, if we take mj = 4.5 GeV and f3 = 6/g 2 = 6.0,6.2,6.4 with a -1 = 
2.3, 2.9, 3.7 GeV respectively, which are typical values for the Monte Carlo sim- 
ulations used in B mesons studies [15], we see that in all cases the logarithmic 
correction is negligible and we obtain the final result: 



f b = 0-83 / 



latt 
B 



(5.6) 
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Note that the prescription of ref. [5] [i.e. x = 1 in eq. (2.11)], would lead to 
much larger values for the renormalization constant Z. As a consequence, one 
would obtain f B ~ (0.97 to 1.07) f B att for typical Monte Carlo slopes of 0.5a -1 to 
0.7a _1 . With the choice x = — 1 [eq. (3.22)] the value of f B would be lowered to 
f B ~ (0.63 to 0.70)/]?". 

Recently the renormalization constants of operators involving heavy quarks 
have been computed [8] by using for the light quarks a nearest- neighbour 0(a) im- 
proved action. In this theory [7] 0(a) corrections are absent from on-shell hadronic 
matrix elements. Using the normalization condition (2.5), the net result for Z 
amounts to replacing the constant 13.59 in (5.5) by 10.08. Then to first order in 
perturbation theory one gets: 

f B =QS7 fatt-imp ^ ( 57 ) 

where f 1 ^ 11 -" 1 ^ i s the Monte Carlo value measured in simulations which employ the 
same improved action as is used in the corresponding perturbative calculations. We 
stress again that in this case the prescription (2.5) [x = in (2.10)] is compulsory 
in order not to loose 0(a) improvement. 

Two numerical calculations of f B with a static 6-quark have been performed in 
the literature [16], [17] at (3 = 6, using for the renormalization coefficient the value 
0.8, essentially in agreement with the estimate (5.6) 4 . A systematic comparison 
with the results on meson decay constants for the case of propagating quarks [3] 
shows a fair agreement with the static points, the latter having however a tendency 
to lie a bit too high. With the EH normalization condition [x = 1 in eq. (2.9)], the 
situation would be much worse, since the static points would be raised by about 
30%. Of course this pattern has to be checked at higher values of (3 5 . 
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4 It is argued in refs. [16] and [17] that the value 0.8 is in essential agreement 
with the results of both [4] and [5]. This statement probably refers to the value 1.22 
found for Z r in [5]. However, as we discussed in section 3, this apparent agreement 
leads in fact to quite different values for f B . 

5 The authors of ref. [17] claim that the static f B decreases sensibly with in- 
creasing (3. 
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Appendix A. 

This appendix is devoted to a sketch of the steps needed to go from eq. (4.2) 
to eq. (4.4). Let us study only l l f a , the analysis for T l £ u being similar. After the 
integrations over po and ko, we get: 



J latt 
1 = 



7 r d 3 p A{$)-Z(z) sr + F d 3 k Zi(fc) 



— (Z<n) T f d k Zl{k) 4 x 

(£)) V 7 (27r) 3 A(fc+£) i-^-Zfi^iffc) 



(2tt) 3 A(£)(^y-Z(£)) _y^ (27r) 3 A( fc +f) ^L_-Z(£)Z l( fc) 

1 l-(^i(fc)) T 



{ 



Zi(fc) 



Zi(fc) l-Z\{k) 
^y-^i(fc) Z(k+£)- Z{z)Z x {k) 



l-(^iW) r ^ m i-(#) 

l-Zi(fc) 



ed-^ } 



(A.l) 



where 



A(fe) = ! + £,- (l-cos(fcj)) 



Zi(fc) = A(fc) - ^(fc) 2 - 1 



Z(fe) = 1 + B(k) - ^/(1 + B(k)) 2 
J2j sin 2 (kj) + (1 - A(k)) 2 



(A.2) 



B(k) 
R 



2A(k) 

Z(k+£) 



and G is the usual Heaviside distribution. 

We are interested in the behaviour of T^ u when r — > oo. It is possible to show 
that the integration over p gives no problem, essentially because it is protected by 
the exponentially damping factor Z(^) T . To illustrate the method we have used, we 
shall not consider the whole expression (A.l), which is rather long, but a simpler 
quantity; this offers however, a degree of mathematical complexity that is similar 
to that of the integrals we encounter in the actual computation, namely: 
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+ 7T 



The problem about eq. (A. 3) is that, since the integral diverges when r — > oo, we 
cannot exchange the limit and the integration. To study the behaviour of J as 
t — > oo, it is convenient to first add and subtract from the integrand its value at 
p = 0; J is then decomposed into a sum J7i + J 2 with: 



' T d 3 k 1 1 l-{Z l {k))'' 



( 2yr ) 3 z(fcy _Zl(fc) zT(k)~ Zl{k) 1 ~ z ^ k ) 



*7i = 



1 1-(Z 1 (^)) T 
^-Z l{ ±) 1-Z l( *) 

Let us call J the integrand in From (A. 2) we see that the quantities Z and Z\ 
are smaller than 1; thus we have the bound: 



i i zfjr -2(^1) i + (!-■?(?>) i i-(ZKj>r 



From (A. 2), after some work, it is possible to show that there exists a set of constants 

1 1 



aj such that 



Z(q) 


Z(q>) 1 
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< ail \q I - I q' 1 1 

(A.6) 



1 _ e -a4k| 
< tt 3 : : 

I q I 

(l-Z(q)) < a 5 \q\ • 

It then follows that the absolute value of J is bounded by an integrable function of 
k and, by Lebesgue's theorem, that we can exchange the limit and the integration 
to get: 

+ f° d 3 k 1 - e~l fc l |fc|-|fc+p|-|p| 
hm J 2 = / ttt-to /i i ;„ 1 3 i r i ; i r , ~i , , , • (A.7) 
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When the actual integral (A.l) is analysed in this way, it turns out that the terms 
that are analogous to J 2 are responsible for the appearance of the constant C\ in 
eq. (4.4). The subtracted terms (analog to J\) give: 



+ F d 3 k Z^k) 2 1 

J (2tt) 3 A(k) i-ZmZ^-Zik) 

x (A.8) 

A(k)-Z(k) l-(Z(k)) T A(k) — Zi(k) 1 - (Z 1 {k)) T 
zfk)-Z(k) l-Z(k) -z^W)-Zi{k) l-Zi(fc) 

which can be rewritten in a simpler form if we introduce an integration over ho'. 

i T d\ e ih ° - 1 - Ej (1 - cos(fcj)) l-e lfc » T 

(2tt) 4 2AiA 2 1 - e ifc ° ' 



-7T 



leading to the first term in the r.h.s. of (4.4). 



Appendix B. 

In this appendix, we shall derive formula (4.7). The zeroth-order graph 1 (a) is 
equal to 3/7r 2 t -3 , and together with eqs. (4.4) and (4.6), we find from our definitions 

c M c t P tree (t) + P s c tat (t) 1 
ptreety + pcty 

„ 4 a ff*k ^-l-E(l-cos(^)) X _^ T 
t-oo 3 f W ( 2 7t) 4 2Ai(A;)A 2 (A;o,k) 1 - e ik ° (B- 1 ) 



°° d 4 " e k e jfcot -l 
(2tt) 4 - ±(fc 2 ) 2 



If we use the identity: 



JL JL , ? sin4( ^ } Ax 

A 2 4Ai + AiA 2 ' A 2 ' ( 2) 
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the lattice integral in eq. (B.l) can be rewritten as: 

+ f d 4 k e ik ° - 1 - E,(l - cos(fej)) 1 - e iko1 

Il = j 



(2tt) 4 2Ai(fc)A 2 (fc ,k) l-e lfc o 

— 7T 

4/ fc 
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+ ~ 1 3 A X A 2 A 2 (0, k) 



1 - e ikoT + F d 3 k 
+ 



l-e* k » J (2tt) 3 2A 2 (0,k) 



When t is large, the factor e lk ° T gives no contribution in the second and the third 
terms in the second equality of eq. (B.3). In the first term, the denominator can 
be exponentiated and one arrives at some combination of Bessel functions. This 
method was discussed in appendix C of ref. [4] and leads to the formula: 



+7T 



MiAr™p( 2h (2) +f - + «) ■ < B - 4 > 

— 7T 

The computation of the continuum integral in eq. (B.l) is straightforward and gives: 



Ic = I 



' (i 4 - e k e ik °> - 1 



oo < 2 -) 4 - (B.5) 
1 1 2 + ln(47r) + 2 In (^-) + lE 



8tH V e v ' V 2 

From eqs. (B.3), (B.4) and (B.5), one finally obtains eq. (4.7) 
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Figure Captions 

Fig. 1: The four contributions of order g 2 to P(t), P s tat(t) and -Pl(t) as denned 
in eqs (1.1), (1-6) and (1.7) respectively (in order to avoid redundancy, we 
use the same graphic representation for the three cases). 

Fig. 2: The graphs contributing to P(t) in an effective-action formulation. 
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